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Abstract 



— . ' We consider ttie Dirac equation in cylindrically symmetric magnetic fields and find 

JjJ I its normal modes as eigenfunctions of a complete set of commuting operators. This 

. set consists of the Dirac operator itself, the ^-components of the linear and the total 
angular momenta, and of one of the possible spin polarization operators. The spin 

' ^ ' structure of the solution is completely fixed independently of the radial distribution of 

the magnetic field which influences only the radial modes. 

' We solve explicitly the radial equations for the uniform magnetic field inside a 

^ , solenoid of a finite radius and consider in detail the scattering of scalar and Dirac 

• particles in this field. For particles with low energy the scattering cross section coin- 

QQ ' cides with the Aharonov-Bohm scattering cross section. We work out the first order 

, corrections to this result caused by the fact that the solenoid radius is finite. At high 

0^ ' energies we obtain the classical result for the scattering cross section. 
O 



PACS numbers: 03.65. Bz, 11.80-m 



1 Introduction 

O . The behavior of relativistic charged particles in external magnetic fields has been the subject 

' of many investigations in QED (see, for example, [|[ D). Most of them has been concerned 

with the synchrotron radiation emitted by charged particles moving in cyclic accelerators 
and storage rings ||l^ . This effect has found wide applications in different fields of physics, 
biology and technology. 

' By its origin the synchrotron radiation is the classical effect. It has been shown that 

quantum corrections became important only for relativistic particles or intense magnetic 
fields [0 0. This happens at high energies E and high field strengths B leading to a 
characteristic product 

B E , M^c^ 

1, with Bo = ^— (1) 



BqMc^ ' ^ eh 

of the order 1. Other QED processes, such as pair production and pair annihilation, have 
no classical counterparts. In these cases too, intensities of the processes become relevant if 
particle energies and magnetic fields fulfil Eq.(|^). Therefore, in QED processes, the external 
magnetic field can not be treated as a perturbation. The relativistic particles must instead 
be described by exact solutions of the Dirac equation in an external magnetic field. This will 
also be done in this paper. For the comprehensive review of exact solutions of relativistic 
wave equations in external fields see Q. 

Much work has been devoted to solutions of the Dirac equation in an uniform electric 
and magnetic fields |7[ ^ and in the electromagnetic plane wave field [pTf. Different QED 
processes have been studied for the corresponding physical situations ||7|, |13|. 

^ e-mail: vdskarzh@sgi.lpi.msk.su 

'^e-mail: Juergen.Audretsch@uni-konstanz.de 
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A new field of research was initiated with the study of the Aharonov-Bohm (AB) effect 
(for the latter see QED processes in the presence of a magnetic string were elabo- 

rated in detail; the differential cross sections of the bremsstrahlung of an electron passing 
by the magnetic string |Q] and the pair production by a single photon in this potential 
have been worked out. In these cases it surprisingly turned out that the cross sections do 
not become small even for low particle energies. This fact can be interpreted as follows: for 
a finite flux the string magnetic field becomes infinitely strong in the limit of a string, so 
that the criterion of Eq.(|l|) is formally fulfilled for all energies. 

To describe realistic situations the magnetic string has to be replaced by a thin solenoid 
which contains an intense but finite magnetic field. It is now very important to know how 
the results obtained in the string limiting case are modified for a solenoid of finite radius. 
We will attack this problem in this and following papers where we are going to investigate 
in detail different QED processes in cylindrical magnetic fields. Below we concentrate on 
the scattering of scalar and Dirac particles. 

In the Sections 2 and 3 we define the normal modes of the Dirac equation in cylindrical 
magnetic field of arbitrary radial dependence as eigenfunctions of a complete set of com- 
muting operators. The set consists of the Dirac operator itself, the z-components of linear 
and total angular momenta, and also, of one of the possible spin polarization operators. 
We fix completely the spin structure of the solution independently of the radial distribution 
of the magnetic field which infiuences only the radial modes. In Sec. 4 we solve explicitly 
the radial equations for the model of a solenoid of finite radius with an uniform magnetic 
field. With the tube radius going to zero this model can be considered as a realistic model 
replacing the AB magnetic string. This approach allows to solve correctly the Dirac equa- 
tion for the AB potential. Then we consider in detail (Sec. 5) the scattering of scalar and 
Dirac particles by the magnetic field of the solenoid of the finite radius. At low energies of 
incident particles the scattering cross section coincides with the AB scattering cross section; 
we derive the first order corrections caused by the finite tube radius. At high energies we 
obtain the classical result for the scattering cross section. 

Throughout we use units such that h = c = 1 . 

2 The Dirac equation in a cylindrically symmetric mag- 
netic field 

The Dirac equation in an external magnetic field reads in cylindrical coordinates (p, ip, z) 

idt-^{x)^ H-^{x), H ^ a,{p, - eAi) + (3M (2) 

where x = {t, p, f,z), e is the charge of the Dirac particle (e > for the positron and e < 
for the electron), ai and P are the known matrices written in the cylindrical coordonates, 
and the kinetic momenta are given by 

TTtp : = - eA^ ^ --d^ ~ eA^p . (3) 

The vector potential for an arbitrary cylindrically symmetric magnetic field of a fixed 
direction (along z-axis) has a nonzero angular component A^{p), and the magnetic field 
reads in terms of the potential 

B,ip)=A'^ip) + - A^ip). (4) 
P 

The dependence of the solutions on z and the azimuthal angle (p can be fixed in de- 
manding that ^!{x) is an eigenfunction of the operators of linear momentum projection 
and total angular momentum projection J3, 

P3*(x) = Pa* (a;), 

J^^ix) = {L3 + 1^3)^ {x)^{~td^ + 1^:3)^ ix)^j3^{x), J3:=l + l-. (5) 
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Here / is the integral part of the half- integral eigenvalue j^. 
We find 



where 

^2 i?2(p) e*^'+'^^ y ' I C4 i?4(p) e'^'+i)'^ 



(6) 



(7) 



i?p = a/p^ + Af = y/p±+I^+M^ is the energy, p_L denotes the radial momentum and 
Cj are the spin polarization coefficients. 

In terms of u— and w— spinors the Dirac equation reads 

aiipi - eAi) u = {Ep + M)v , 

{p, -eA{)v = {Ep - M)u . (8) 

From Eqs.(|^) and we find equations for the radial functions 

\ P 

-iCiR[ +i(- 
\P 

+-1 



P 



eA^^ 


1 C2i?2 


+ P3C1R1 = 


(Ep 


+ M)C3R3 , 


eA^^ 




- P3C2R2 = 


{Ep 


+ M)C4,R4: , 


eA^^ 


j C4R4 


+ P3C3R3 = 


{Ep 


- M)CiRi , 


eA^ 


I C3R3 


- P3C4R4: = 


{Ep 


- M)C2i?2 . 



-iC3i?3 +i(^^-eA^] C3R3 - P3C4R4 = {Ep - M)C2i?2 . (9) 
The spin coefficients can be constrained by the normalization condition 

i=4 



Eicd' = i- (10) 



3 Spin polarization states 

In the case of cylindrical magnetic fields the spin polarization coefficients C'i can be fixed 
by one of two spin operators St or S3 which commutes with operators H,P3 and J3 but 
do not commute with each other. The hclicity operator, 

S. := 5:*^, (11) 
P 

describes longitudinal polarization (the projection of the spin onto the velocity direction of 
the Dirac particle). The z-component S3 of the operator S := f3 T, + j{p~ eA)/M , 

^3-/^^3+7|, 7=(; ;) (12) 

defines a transverse polarization (along the direction of the magnetic field) for nonrelativistic 
motion, or for the motion in the plane perpendicular to the magnetic field. The operators 
St and 5*3 are not independent. However, to make things easy for the reader we will discuss 
them separately. 

3.1 The helicity states 

For the solution of the Dirac equation (||) which is the eigenstate of the helicity operator 

5t^'(a;) = s^'(x), s = ±1, (13) 
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the Eqs.dq) and (O) connect the u— and w— spinors according to 



v = s ^ , ^ =u. 14 

It means that ~ Ri, = R2 and that the coefficients 0^,0^ are coupled by Eq. ([l^ ) 
to Ci , C2 correspondingly. Choosing Ci and C2 to be 



^£:p + M^p + sp3 „ isy/Ep + M^p- sp3 

and determining coefScients C3 and C4 from Eq.(|lJ), 

^3 — , — = , 04 = , — -= , (iol 

we obtain from Eq.(^ the following set of the equations for the independent radial compo- 
nents i?i,2, 

^2 + ("~~" ^ ^^v^ ^2 = P±Ri , 

-R[ ^1 = P^R2 . (17) 

The solution of the Dirac equation with the quantum numbers Ej p^jj^, s has then the form 
of Eqs.(|), (0) with R3 = Rj^,Ri = R^ and the coefficients (|l|), (|l|). 

3.2 The transverse polarization states 

On the other hand, for the eigenstate of the operator 53 (p^), 



S3^{x)=S3^ix), S3^S\/1 + ^, s = ±l, (18) 



we find from Eqs.(|) and (y) that 

^^1 = 777 m^i' '"2 = 777 T\'^2. (19) 

M[S3 + 1) M(s3 - 1) 

This means that R3 = Ri, R^ = R2 and that the coefficients C3 , C4 are defined by (|l|) 
from Ci , C2 , correspondingly. Fixing the coefficients Ci and C2 by the expressions 



^^ = \ -^ 7S^= — ' ^ ie(S3P3)\ — == — (20) 



2s3 ^/2e: 'V 2s3 



where e{x) is the sign of x and determining the coefficients the C3 and 6*4 from Eq.(|l9|), 



we obtain from Eq.(^ again the radial equations (p^). The solution of the Dirac equation 
with the quantum numbers E, Ps^.ls: S3 has therefore the form Eqs.(|^), (^) with i?3 = 
i?i , i?4 = i?2 and the coefficients (|20| ) , (^) . 

We have shown how the spin polarization operators (|ll|) or (^2|) fix the spin structure 
of the solutions of the Dirac equation in arbitrary cylindrical magnetic field. Note that 
solutions for different configurations of the magnetic field differ from each other only by 
their radial functions. 
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3.3 The charge conjugate states 

The solutions ^ of the Dirac equation (||) represent electron (e < 0) and positron (e > 0) 
wave functions of positive energy, '^e,p{x). The complete set of solutions of the Dirac 
equation (||) includes the positive and negative energy electron (or positron) states. Instead 
of negative energy electron (or positron) states one can use as well the charge conjugated 
positron (electron) states which can be obtained by the charge conjugation operation, 

*(x) *e(x) := C^transp(x), C = ^2 . (22) 

For the helicity and transverse states we obtain 



^,{x) = ^ e^E,'-^^^. i>.{pM, Mp,v) - f "0 



where 



iCl RUp) e-*('+i)'^ 
~iC^ Rlip) e-'^f 



iCl Rlip) e 



(23) 



(24) 



with the coefficients (p^, (16) and (pO|), (pl|), correspondingly. 

With reference to these states, the electron-positron field operator reads 

4^{x)^Y.^^,ix)a,+%ix)b]] (25) 
j 

where Uj and bj are the annihilation operators for electrons and positrons with quantum 
numbers j := {p±,P3, 1, s). This field operator obeys Eq.(||) with e < 0. Of course, we could 
use the charged conjugated field operator which obeys this equation with e > 0. 



3.4 The radial equations 



The set of the radial equations (jl^) is equivalent to the following second order equation for 
the first component 



R'l + -R[ 
P 



- - eA. 



i?i -0; 



the second component is then defined by the equation 



i?2 



1 



eA^ Ri 



(26) 



(27) 



Alternatively one could start with the second order equation for R2 and define the first 
component i?i by the equation (|lj). Both ways are completely equivalent. We would like 
to stress that because of the first order constraints ( p^ ) between the components i?i and R2 , 
the Dirac equation differs from the Klein-Gordon equation for relativistic scalar particles. 
The radial functions can be normalized by the condition 



/ 



pdp Rl 2{p'±P) RiAp^p) 



s{p± -p1) 
\/p±p'± 



The Eqs. (|26|), ( p7|) together with the normalization condition (^ 
tions Ri and R2 up to a common phase factor. 



(28) 



define the radial func- 
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The nonrelativistic limit: The positron and electron radial and angular functions 
conserve their forms in the nonrelativistic approximation {Ep = M + £p, £p ^ M) with 
nonzero spin coefficients 



Ci = = — , G2 = 1= — , tor liehcity states 

V2p 



Vp 


+ SP3 


V 






+ 1 



C\ = ^ "r— " ; ^2 = — transverse states (29) 
'2s y2s 



4 Solenoid of a finite radius and the magnetic string 

Our basic aim is to study the AB scattering of scalar and Dirac particles for experimentally 
realistic situations. We consider therefore a solenoid of finite radius po with the uniform 
magnetic field inside. This model contains the AB magnetic string in the limit of zero radius 
and constant magnetic flux. This limit can alternativly be based on a different model |^ 
for which the magnetic field is concentrated on the surface of a cylinder of radius po- It is 
obvious that this second model being more simple for the investigation of zero radius limit 
does not describe a realistic experimental set up. 

4.1 The uniform magnetic field inside a solenoid of finite radius 

The uniform magnetic field (0, 0, B) restricted to the interior of a solenoid of finite radius 
Po can be described by the vector potential 

\ e$/27rp \ (j)/p li p>po ^ ' 

where $ — ttPqB := /3$o ~ ^fp^o is the magnetic flux in the solenoid, $0 := 27r/|e| is the 
magnetic flux quantum and e is the sign of the charge (e > for positrons and £ < for 
electrons). We take /3 > so that the constant magnetic field inside the solenoid 

B = ^, P<Po (31) 
\e\PQ 

points in the positive z direction. To obtain the uniform magnetic field in the whole space 
one may perform the limiting procedure f3 \e\BpQ/2 together with pa 00. 

Let us consider first the internal solution. The general radial equations (^6|) and (|21 
read in this case 



R'l H — R'l 
P 



i?2 



2Pj -^—-PJL 

P Po J Po 



1 

P± 



Ri+{-~^p]R 

P Po 
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i?i=0, (32) 
(33) 



Introducing the new variable x ~ P{p / Po) and the new functions ^1,2(2;) according to 

Ri{x) = e"* Yi{x) , i?2(x) = x-^ e-t Y2{x) (34) 



we obtain from Eq.(|32|) the confluent hypergeometric equation |g] 

xYl' + {Bi - x)Yl - AiYi ^ Q , (35) 

Ai = t£i±l^-^, Bi = \l\ + l. 
2 4/3 ' ' ' 



Then from Eqs.(|33|) - ( p5[ ) we find the following confluent hypergeometric equation for the 
function Y2{x) 

xYi' + {B2 - x)Y^ - A2Y2 = , (36) 

A2 = M_£(i±l)±l±£_d|l, B2^\l + l\ + l. 
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Solutions of Eqs.(^5|) and (^6|) which are regular at p = are the confluent hypergeometric 
functions ^{Ai, Bi; x) and $(^2, -B2; a^), correspondingly, so that we obtain 



Riix) = ci X- e-^ ^{Ai.Bi ; a;) , 



B'^^x) 



c; a: 2 e 



$(A2,B2;a:) 



(37) 



The coefficients q and c; are connected with each other by the Eq.(|3£ 

The external solution obeys the radial equations ( ^6|) and (^) at p > po- 



1 



R'l H — R'l 



i?2 = — 



Ri+pI Ri^O, 



I - 



-Ri 



(38) 



Solving these equations by Bessel functions of the first kind with positive and negative order 
and using Eq.(p7|) we find the external radial components 



where 



R^ip) 
RTip) 

1^ 



ai Juip^p) + bi J-^{p_i_p) , 
e/_0 [ai Ji,{p^p) - bi J_o{p^p)] 



1 

-1 



if / > 
if / < 



(39) 



(40) 



The coefficients a/ and 6/ can be obtained from the matching conditions at the surface 
p = po or (cc = /3) of the solenoid for the internal and external solutions. Since radial 
components are connected with each other by Eq. (^3|) one can use any couple of the matching 
conditions for components R\ , R2 and their first derivatives R!-^ , i?2 . Choosing the conditions 
for R\ and R!-^ we obtain 



Ci/?^' e 2 $(Ai,Bi;/3) ^ai J^{pA_po) + bi J-^{p^po) , 



(41) 



P±Po [ 2/3 $(Ai,Bi;/3) 



c;/?'^ e 2 $(Ai,Bi;/3) = aiJlip^po) + biJL^ip^po)- 



Eliminating the coefficient c/ from these equations we find the matching conditions in the 
form 



2f3 
P±Pn 

w 

P±Po 



\l\-(3 , $'(Ai,i?i;/3) 



2/3 $(Ai,Bi;/3) 
\1\~I3 , ^i$(Ai + l,Bi + l;/3) 



(42) 



2/3 



Bi $(Ai,Bi;/3) 



ai JliP±Po) + bi JL^ip±Po) 
ai Ju{p±Po) + bi J^u{p±Po) 



which fixes all coefficients up to a normalization constant. The same results can be obtained 
for instance from the continuity conditions for Ri and R2 at p — pQ. 



4.2 The zero radius limit (magnetic string) 

The vector potential for the infinitely thin, infinitely long straight magnetic string lying 
along z-axis 

e$ e$ £/3 (j) 



eA^ 



< p < 00 



(43) 



27rp $oP P P 

can be obtained from the vector potential of the solenoid ( pO[ ) for po — > keeping the 
magnetic flux constant. It is singular on z-axis and produces the singular magnetic field 



B = ^S{p) 
\e\p 



(44) 



7 



which is concentrated on the z-axis and points to the positive z direction. We separate the 
integral number N from the flux parameter P, P = N + 6, N > 0, < 5 < 1 since it is 
the fractional part S of the dimensionless magnetic flux which produces all physical effects. 
Its integral part N will appear as a phase factor exp{ieN(p) in the solutions to the Dirac 
equation. 

In the case of the AB potential the radial solutions are given by Eqs.(|39|), but have an- 
other domain of definition. We note that the component i?i contains modes with coefficients 
bi which are singular at p = for all I. The component R2 contains singular modes with 
coefficients qn (e > 0), a_jv_i (e < 0) and &;^jv (e > 0), 6i^_jv-i (e < 0). The appearance 
of the singular modes is usually forbidden by the normalization condition (|2^ ) for radial 
functions. But for unbounded potentials, and for the AB potential ( ^3|) in particular, this 
condition is less restrictive. The normalization condition ( p8|) requires that the integrals of 
the type 

/•oo 

I^. = / 4{P±P)pdp (45) 
Jo 

were convergent at small p what takes place at /i > —1. For the positron solution (e > 0) 
the normalization condition (|45| ) for Ri allows nonzero coefficients ai, bjy and 6jv+i while 
this condition for R2 allows nonzero coefficients a;, 5jv and bj^^i. For the Dirac equation 
both components Ri and R2 must satisfy the normalization condition (^5|). It means that 
only the coefficient b^ can be nonzero presenting one singular mode with I = N in the 
positron solution of the Dirac equation (while two singular modes are allowed for scalar 
wave equations). For the electron solution (e < 0) the normalization condition (^5|) for Ri 
allows nonzero coefficients ai, 5_jv and b-N-i while the condition for R2 allows nonzero 
coefficients ai, b-N-i and 6_jv-2- Therefore only one singular mode with I = —N — 1 can 
be present in the electron solution of the Dirac equation. 

The problem of the singular modes is related to the fact that the Dirac operator as 
well as any Hamilton operator for charged particles is not selfadjoint in the presence of 
the AB potential. This would cause many problems for the unitary evolution of quantum 
systems unless selfadjoint extensions of these Hamilton operators exist. The selfadjoint 
extension procedure applied to this problem gives results which can be obtained by the 
direct calculation of the normalization integrals (for a detailed discussion see [^). However, 
this procedure does not fix the extension parameters which determines the behavior of the 
wave function at the origin. This situation is not satisfactory from the physical point of 
view. It can be overcome in turning to better defined models in particular, to the model 
of a solenoid of a finite radius with the uniform magnetic field [Q. 

From Eq.(|4^) we obtain for po ^ 

Bi ^[Ai,Bi;P) a + 6 

(47) 



where 



= ( p^Po Y" r(-t^ + 1) 



and 



, u-\l\+P-2p {A,/B,) (MAi + 1,B, + I; p)/'P{A,, B,; P)) 

a/ ^"j^+|Z|-/? + 2/3(Ai/Si) ($(Ai + l,Si + l;/3)/<i>(Ai,i3i;/3))- ^ > 

The asymptotic behavior of S,i for po — » depends on the behavior of the denominator in 
Eq. (^) . Since > we have bi ^ Pq" ^ for po ^ unless the following equality is valid, 

■in. D:^ to (.+ |,|-, + 2,f it|±i.f±i:«) .0, (49) 

PO— >oo po^oo Y Bi <P(/ii, _Di; pj / 

This can happen only if 

iy+\l\-l3^0 and lim Ai ^ \l\ - el + 1 - e ^ . (50) 

po^oo 
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In this case the parameter S.i goes to zero, 

6 ~ Po"^^'' ^0 at po ^ oo , (51) 

unless 

ly ^ \l ~ < I. (52) 

The Eqs.(|50|) and ( [5^ ) are fulfilled for the positron (e > 0) singular mode I = N only. It 
means that aN = 0, 6jv 7^ 0, and that accordingly the first positron component with I = N 
is singular and the second component is regular at p = : 

Ri{p) = J-s{p±p) , R2{p) = Ji-s{p±p) ■ (53) 

In this case the interaction between positron magnetic moment and the string magnetic 
field is attractive for the first (upper) component of the wave function and it is repulsive 
for the second (down) component. For the electron (e < 0) solution all coefficients bi ~ 0. 
It means that the singular electron mode / = — A'^ — 1 has the regular first component and 
the singular second component 

Riip) = Ji-s{p±p) , R2{p) ^ -J-s{p^p)- (54) 

Introducing new notations we can rewrite the positron and electron solutions of the 
Dirac equation in the presence of the magnetic string in the following final forms which are 
valid both for regular and singular modes, 

R,{p) = e^tl'-^^l J,,(p^p), 

R2{p) = e^^l'-^^l e,_,Ar J,,(pxp) = -ie*5l'+i-^^l J,,(pip) (55) 

with 

vi := ei^^N {I - <f>) , 1^2 ■= f^i-eN {I + 1- cj)), (56) 
_ / 1 if ; > eiV 
''-^^ I -1 if ; < eiV ■ 

These radial solutions satisfy the normalization conditions (^). 

Redefinition of I : The matrix elements of the QED processes (bremsstrahlung 
of electrons and positrons, pair production and annihilation) in the presence of the 
AB potential contain integrals over the products ^'(a;)7^A^(a;)^'(a;), ^'c(a;)7^A^(j:)\E'c(x), 
^'(a:)7p^^(a;)\E'c(a;), \E'c(a:)7^A^(a;)\E'(a;), correspondingly, where '^{x) and "^dx) are the 
electron and charge conjugate positron functions. The integer number N disappears from 
all these matrix elements after the redefinition of the angular quantum number I, I l+eN. 
This means that the matrix elements of the QED processes are independent of the integer 
part of the magnetic flux in units of the magnetic quantum. One can foresee this fact 
beforehand since the Dirac equation (||) with [3 = N + 5 can be transformed to the Dirac 
equation with j3 = 5 hy means of the gauge transformation 

*(a;) ^'{x) ^ eyiY>{-ieNif)^{x) , 

sfj sip — N) 

eAJx) :— — — > eA'Jx) — eAi(x) — i exp(zeiVi^)Vi exp(— ieiVi^) = . 

P P 

All observable quantities are conserved under this gauge transformation. 



5 Scattering of scalar and Dirac particles by a solenoid 
of finite radius 

The expressions obtained above present the partial positron and electron wave functions in 
terms of cylindrical modes. These states do not describe outgoing particles with definite 
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linear momenta at infinity. In order to calculate the cross section of QED processes we 
need the positron and electron scattering wave functions. In external fields there exist 
two independent exact solutions of the Dirac equation, '^^^\p; x), which behave at large 
distances like a plane wave (propagating in the direction p given by Px = p± cos (pp, Py — 
p± sin(pp, pz) plus an outgoing or ingoing cylindrical wave, correspondingly. Because of the 
damping of the cylindrical waves at large distances we may use these superpositions instead 
of plane waves. To evaluate correctly the matrix elements in the presence of external fields 
we have to take wave functions for ingoing (outgoing) particles which contain outgoing 
(ingoing) cylindrical waves. We turn first to a discussion of scalar particles. 



5.1 Low energy scattering of scalar particles 



The wave function for scalar particles can be found as the eigenfunction of the complete set 
of the commuting operators (o): 



Its radial modes i?;(p) obey the radial Klein-Gordon equation. 



(57) 



(58) 



which is similar to the Eq.(p6|) but does not contain the term e(A^ + A^j p) which describe 
the spin-magnetic field interaction. 

For realistic models the internal radial solutions are some regular functions c; i?j"*(p). 
For model ( ^ ) of a solenoid with constant magnetic field this is 



Q x 's e 2 S; x), x = ^ , 



(59) 



I + 1 PjPo 
2 4/3 



Bi 



The external solutions are similar to the first line of Eqs.(|39|). It is convenient to rewrite 
them as follows 



i?/(p) = (a/ + 6, e "'^) J^(j3ip)+i6/sin7ri/i/W(pip), 



(60) 



where h'^^ (x) are the Hankel function. The corresponding scattering wave functions can 
be obtained by the superposition of these cylindrical modes. 

Since the plane wave term of the scattering wave function at p — > cx) is defined only by 
the external solution, the scattering wave function can be written as follows 



1 



^{p-x)^ — e-''^-'^'^^^W.P.V). 



(61) 



Jltfi _ 



J\l^eS\iP±P) - Plifp) hI-1}^Jpj_p) 



jK'-p-'-Pp) 



with the coefficients 



The wave function (^) behaves asymptotically as follows 



V'(p ; p, if) ~ e'^^^'^-'^"-'^) X 



(62) 



(63) 
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where 



V^vrpj_ ^ L J 

is the scattering amphtude. We note that the distorting factor e'^^'-'^"'^''"'^^ appears in 
Eqs.(|6^) and (64) because the vector potential (30) decreases slowly at infinity. 

The wave function (^TJ) and the scattering amplitude (Hh contain an arbitrary coefficient 
Hii^Pp) connected with the coefficients a; and 6; of Eg. (|60[) . It describes how the finite 
radius pq of the solenoid influences the scattering cross section and can be found from the 
corresponding matching condition: 

Mp±Po) {lnMp±Po)y - i\nR]-\p^po)y 

pLi-eNKVp) = — m m ■ • (65) 

hI'\p^Po) {In Hi'\p^po)y - {\nR\-'{p^po)y 

The partial wave decomposition is the effective method for the investigation of low energy 
scattering. In this case the de Broglie wavelength of scattered particles is large compared 
with the size of the target (pa), and the lower angular momenta contributes mainly to the 
scattering cross section. One can see from Eq.(|6^) that pi ^ at po ^ 0. For any realistic 
model, including the model (po[), the internal solution is regular at po = 0. Then we have 



^ L._.. . hm (ln^.b.Po))--(lni^rb.Po))- ^ ^^^^^ 

Ri[P±Po)) P±Po Po-^" {\nHl)'{p^po)y - {lnR\'^\pj_pQ)y 

and the second term in Eq.(|^) disappears completely at po 0, 



/X;-eAf(v3p) ~ ZSmTTt/ r(t/ + 1) ^'"''^ V~2~/ ^ ^' ^^"^^ 

with 

,^ ^ \l\ + P - 2(37ji ^'{AuB,; P) 

Accordingly we have obtained the following result: At low energies, pxpo 0, the scattering 
wave function (^Tj) coincides with the AB wave function which behaves asymptotically as 
follows 



(69) 



in whole space outside a narrow region around the forward direction. The corresponding 
AB scattering amplitude reads 

1 .„ . v-fp sin7r(5 

fABiip - ^p) = -e -= e-^ (70) 

and the differential cross section of the AB scattering m is equal to 



da p± 2 1 sm nS 

dip p 2'Kp sixT ^ 



To find the first correction to the AB scattering amplitude at small po one need to take 
a minimal value of i^. It is (5 at < (5 < 1/2, I ^ eN and l-(5atl/2<(5<l, ; = e(A^+ 1). 
We obtain 

( \ ■ ■ . r(i-^) /piPo\2^ 

p^[ipp) = ismird Yi^i + S) ) ' 

t ^ ■ ■ . r(J) /Pipofd-*) 

Aie((pp) = zsuiTTd L, ) (72) 
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with 



N+S'U(N+l) 



{N + 5)i^,N ■ - 1 
Then the scattering wave function reads as follows 



(73) 



= (7V+l)(ip-V=p-7r) 



(74) 



where ipAsiP', P, ^) is the AB scattering wave function with the asymptotic behavior (69) 
Now the additional term changes the modes with I — eN, e{N + 1) of the scattering amph 
tude fAsiv - Vp) given by Eq.(|7^) 



Sfif ~ fp) 



fAB{(P - (Pp) + 5f{ip - ipp) 

2 



^/2^^pZ 

In this approximation the scattering cross section is equal to: 
d(T P±,„/ m9 1 sin^ 7r5 , 

- y'/(---^)'^-2^^i^?s(i + ^)' 

A = 4sin^^sin(^^^+.£<5) 



2 ' 2 

r(i + ^) °V 2 / ' r(2-(5) 



(76) 



whereby A = for po = 0. 



5.2 High energy scattering of scalar particles 

The partial wave method becomes ineffective when applied to the scattering of high energy 
particles. The smaller the de Broglie wavelength of scattered particles, the higher angular 
momenta contribute to the scattering cross section. At p^_po oo this corresponds to 
I ^ oo. In this case the scattering process can be expected to be the quasiclassical one and 
the WKB method ||ll| can be used to obtain the scattering cross section. 
The quasiclassical radial wave function for the model (pQ) is equal to 



WKB 



= const " (p) exp 



Then the phase shift in the WKB approximation is given by 



^WKB 



dp' VF]{p)-pi_ 



P±Pn 



(77) 



(78) 



where we introduced the radius of the orbit in the uniform magnetic field R — p±/eB. 
Mmin is the classical turning point of the motion inside the solenoid where the radial kinetic 
energy Fi(p) goes to zero. Fi{p) contains pQ via A^, and in the zero radius limit Eq.(|7^) 
agrees with the corresponding expression for AB phase shift. 

The classical deflection angle for charged particle moving in the magnetic field of the 
solenoid of finite radius ( ^o| ) can be obtained from the corresponding radial Hamilton-Jacobi 
equation. The WKB method gives the same result. Differentiating the phase shift (fz^ ) with 
respect to I we obtain the classical deflection angle 



WKB 
I 



= TT- 2 



dp 



dp 



P±Po 
2R 



2R 



(79) 
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Introducing the impact parameter 



2 

2R 



and calculating the integral (um we obtain 



tan^ = #HZ. 
2 R + a 



(80) 



(81) 



One can see from Eq.(^l]) that one or two values of the deflection angle correspond to 
one and the same value of the impact parameter a. In this case classically unobtainable 
interference effects arise in the framework of the WKB approximation. Since the difference 
between the two impact parameters that lead to the same deflection angle is of macroscopic 
size in comparing with the particle wave length these interference effects are unobservable 
and can be neglected. 

Then we obtain that the differential scattering cross section of scalar particles by the 
magnetic tube of the finite radius pg in the WKB approximation is equal to the classical 
cross section 



dip 



da 

dip 



in 

r sm ^ 



\/Po^R'^ sin 



2 !£_ 

2 



sinf ^ 



aX R < pa 



aX R > Pa 



(82) 



Integrating over p we find that the total scattering cross section of high energy incident 
particles is equal to cr = 2po- This means that the scattering cross section decreases with 
the energy of incident particles approaching its geometrical value. 



5.3 Low energy scattering of Dirac particles 



In this section we consider the scattering of Dirac particles of low energies by the solenoid 
of the finite radius po- As the initial approximation at po — > we prefer to use the radial 
function in the form of Eqs. (|55|), (|5^ ) instead of Eqs.(p9|). With this choice we avoid the 
problem connected with the specific behavior of the parameter ^/ of Eq. (^ ) at po ^ with 
I — N ioT the positron states and I — —N — 1 for the electron states. 

Then, for the solenoid of the finite radius, the external positron and electron radial 
solutions read as follows 



Riip) 

R2{P) = 



(83) 



with the orders of the Bessel functions given by Eq. (^6|) . The coefficients p/ arc defined by 
the corresponding matching conditions and decrease rapidly for po — » 0. With these radial 
functions the cylindrical modes for the Dirac equation (|^) are given by Eqs.(|l) and (0) with 
i?3 = i?i,i?4 = i?2 and the spin coefficients (|T5|), (^ or (pO|), ( pi] ) correspondingly. 
The scattering wave functions for the Dirac equation (0), 



'^I^d{p;p,p) 



u{p;p, pY 
v{p;p,ip) 



are defined by the series 



sN 



CiRi{p)e 
C2 R2{p) e 



dip 



sN 



C^Ri{p)e 
Ci R2{p) e 



Ci ipiiP^Pi v>) 
C3tpi(,p;p,'p) 

-id e*"^" tp2{p;p,p) , 



(84) 

(85) 
(86) 
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with the coefficients 

These coefficients differ from the coefficients h-eNi'^p) ( p^ of the scattering wave function 
for the scalar particles (|6^) only at Z = iV for the positron mode and at I — —N — 1 for the 
electron mode. 

Using Eqs.(|83|) we find the components of the Dirac scattering wave function for the 
positron solution 

■ip2{p;p,(p) = ■ilj{p;p,(p) + Silj2{p;p,(p) (88) 
and for the electron solution 

= ip{p;p,ip) + 5il}i{p;p,ip), (89) 
Mp: p, if) = V(P ; P, ^) + e-^^('^-'^-— ) i smnS e'^' H^^ (p^p) + SMp; P, ^) 

Here the function •i/j{p; p, f ) is the scattering wave function ( [6l| ) for scalar particles, and 
the terms of the first approximation are equal to 

SMp;P,^) = -e^^('^-'^''-^)^e'^l'l e^i^'V^<(U)(Pi/')e''^'^"''^^ (90) 

I 

for the positron solution, and 

6Mp;P,^) = -e-*^('^-'^''-")^e'tl'l e*f(l'l-l'+*l) i/W^,|(p^p) e''('^-'^''\ (91) 

I 



for the electron solution. These terms decrease rapidly for po — > 0. However Eqs.(|8|) and 
( jsll ) contain the additional terms which appear in the AB scattering wave functions for 
the Dirac particles for po and which are absent for scalar particles. They arise due 
to a distortion of the positron radial mode with I = N and the electron radial mode with 
I = —1 — N caused by the interaction the magnetic momenta of the Dirac particles with the 
string magnetic fields. This interaction does not influence other modes which go to zero on 
the string. 

At p±p — > oo, and taking into account Eqs.(|69|), we obtain the asymptotic behavior of 
the Dirac scattering wave function (p4) - (p9h. 



ip{p;p, (p) 



u{p) 
v{p) 



(92) 



VP 

where the scattering amphtude matrix ^'(9:' — (fp) G{(p — (pp) is defined by the expressions 



and 



Sf{^ - Vp) = \ e"'^ y\ e^""'^ p.1 e^'(^-^-) (93) 

V 7rp_L ^ 
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with 

D{V-^p)=\ Q ^. ^--^p ) = cos — 1(73 sm — 7, — • (95) 



The coefficients /i; in Eq.(93) are defined by the matching conditions and go to zero at 
Here the plane wave bispinor 



describes the ingoing particle with the linear momentum p. 

The differential cross section of the scattering for the Dirac particles is equal to 

where p j and jp are flux densities of the scattered cylindrical and the incident plane waves. 
Since G'^{(p — (pp) cxpif ) G{(p — ipp) — ap{ipp), and the relations 



^\p)ap~i^p)m = i>Hp}(—Mvp) + -c^^m = 4- (98) 

VP P / ^p 

are valid both for the helicity and the transverse states, we find that the differential cross 
section of the scattering of positrons and electrons by the solenoid of the finite radius 

is independent of the spin polarization quantum number s. It coincides with the differential 
cross section for scalar particles. This means that spin effects do not appear in the scattering 
of spin particles in given helicity or transverse polarization states which are conserved in 
this process. 

For particles polarized along an arbitrary vector n and described by the density matrix 

TZ{n) = ^{l + a-n), (100) 

the differential cross section can be calculated by the formula 

^ = ^ \F{ip - ^p)\^Tt n{n') D{ip - ipp) n{n) D^ip - ipp) (101) 
dip p 

where n' is the direction of the polarization vector after the scattering. 

Calculating the trace on the right hand side of Eq.(lOl) we find the differential cross 
section for particle beams of the given polarization 

^ = — \F{tp-(pp)\'^ ■^{{l+nn') + {nn'-n3n'^) cos{tp-ipp) + [fixfi']3 sin(<y9-(^p)}. (102) 

For the case of the same polarizations of the scattered and the incident beams, i.e. at 
n = fi' , we obtain 

I - ^ \Fi^ - ^pr (cos^ ^+nl sin^ ^) . (103) 
One can see that the spin effect occurs for polarized particles beams. 
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6 Conclusion 



We studied in detail the Dirac equation in cylindrical magnetic fields of a fixed direction 
and found eigenfunctions of a complete set of commuting operators which consists of the 
Dirac operator itself, the ^-components of linear and total angular momenta, and of one of 
two possible spin polarization operators. The spin structure of the solutions is independent 
on the radial distribution of the magnetic field and is fixed by one of the spin polarization 
operators. The magnetic field influences only the radial modes. We solved explicitly the 
radial equations for the model where the magnetic field is distributed uniformly inside a 
solenoid of a finite radius. From this realistic model we obtained correct solutions of the 
Dirac equation for the AB potential (magnetic string) in the zero radius limit. Also we 
considered carefully the scattering of scalar and Dirac particles by the magnetic tube of 
a finite radius. For low energy particles the scattering cross section coincides with the 
AB scattering cross section. In this case the de Broglie wavelength of particles is bigger 
compared with the tube radius and the influence of enclosed magnetic flux leads to the large 
cross section. This mechanism of interaction of quantum particles with the string magnetic 
field is more effective than the local interaction with the magnetic field. We calculated 
the first order corrections connected with a finite radius of the tube. When the energy 
of incident particles increases the influence of the enclosed magnetic flux weakens and the 
cross section decreases. At high energies we obtained therefore the classical result for the 
scattering cross section. 

The exact solutions derived above will serve as well as basis for future calculations. We 
plan to study how the flnite size of the solenoid influences the differential cross sections of 
the bremsstrahlung emitted by an electron and the pair production by a single photon at 
different energies of the incoming particles. 
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